Superconductivity in molecular solids with Jahn- Teller phonons 
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We analyze fuUeride superconductivity at experimental doping levels, treating the electron- 
electron and electron-phonon interactions on an equal footing, and establish the existence of novel 
physics which helps explain the unusually high superconducting transition temperatures in these 
systems. The Jahn- Teller phonons create a local (intramolecular) pairing that is surprisingly resis- 
tant to the Coulomb repulsion, despite the weakness of retardation in these low-bandwidth systems. 
The requirement for coherence throughout the solid to establish superconductivity then yields a very 
strong doping dependence to Tc, one consistent with experiment and much stronger than expected 
from standard Eliashberg theory. 



The discovery of superconductivity in alkali-doped 
Ceo, persisting up to unexpectedly high temperatures 
(Tc = 33 K [ll or = 40 K |]), raises interesting ques- 
tions about superconductivity in low-bandwidth molec- 
ular solids. Superconductivity arises from an effective 
attractive interaction between the electrons. In conven- 
tional superconductors a net attractive interaction sur- 
vives, in spite of the strong Coulomb repulsion, thanks 
to retardation effects However, retardation is small 
for the fuUerides Q,^, since the molecular vibration fre- 
quencies are comparable to the bare electron bandwidth. 
We show that the combination of molecular solid char- 
acter and coupling to Jahn- Teller phonons produces a 
local pairing, important for superconductivity, which is 
not strongly suppressed by the Coulomb repulsion. In ad- 
dition, the transition temperature depends anomalously 
strongly on the doping level. The superconducting mech- 
anism in fuUerides therefore differs in important ways 
from that of conventional superconductors. 

Conventional superconductors are studied in the 
Migdal-Eliashberg theory, assuming a band width W 
much larger than a typical phonon frequency LUph- For 
the fuUerides, cuph ^ W, so the Eliashberg theory is of 
questionable accuracy. This failure of Eliashberg theory 
is typically thought to lower Tc |^ (although the oppo- 
site has also been argued Q]). Metallic fuUerides have 
very large, nonsaturating resistivities in the normal state 
[H,D suggesting "bad metal" behavior [|o| which is also 
expected to reduce Tc jl^. However, we find that Tc in 
the fuUerides is not generally lower than expected from 
Eliashberg theory due to an unusual cancellation of coun- 
tervailing effects. The violation of Eliashberg theory as- 
serts itself explicitly in a very strong doping dependence 
of T,. 

In AaCgo (A= K, Rb), the three-fold degenerate tiu 
level is partly occupied and couples strongly to eight Hg 
intramolecular Jahn- Teller phonons. We capture the es- 
sential physics using a model with one ii„ level and one 
Hg mode per molecule, with a dimensionless electron- 
phonon coupling strength A. We also include the hop- 
ping between the molecules and the Coulomb repulsion 
U between two electrons on the same molecule llTTll . The 



model explicitly includes Jahn- Teller coupling and places 
no implicit restrictions on the ratio ujph/W or the value 
of A. We refer to this model as the T x h problem. To 
reveal the novel effects of Jahn- Teller character, we com- 
pare this model to a nondegenerate (a) or two-fold degen- 
erate (e) level interacting with a non-Jahn- Teller Ag or 
two- fold degenerate {Eg) phonon, i.e., T x a, E x e and 
A X a problems, respectively. Typical parameters are 
A ~ 0.5 - 1, LUph/W ~ 0.1 - 0.25 and U/W ~ 1.5 - 2.5 
II- 

We circumvent the limitations of Eliashberg theory by 
using the dynamical mean- field theory (DMFT) (l|] with 
a non-perturbative Quantum Monte-Carlo (QMC) tech- 
nique iQ. The electron self-energy is assumed to be q- 
independent, allowing a mapping of the lattice problem 
onto an effective impurity problem. We study supercon- 
ductivity by applying a perturbation creating electron 
pairs and calculating the corresponding response func- 
tion, i.e. the q = pairing susceptibility x- A divergence 
of X below a temperature Tc signals the onset of super- 
conductivity jl^. We write 



x = (1-xor) 



Xo, 



(1) 



where xo is obtained from products of two fully dressed 
electron Green functions describing the propagation of 
two electrons (holes) , which do not interact with each 
other. Eqn. (P then defines the effective interaction F. 
We define a local (intramolecular) pairing susceptibility 



(2) 



where (...) denotes a thermal average and m labels the 
tiu levels on one molecule. Then 
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X^"" and Xo°'' can be calculated within DMFT; this de- 
fines the local interaction F'°'^. F w F'""^ should be a 
rather good approximation, since the interaction is dom- 
inated by intramolecular phonons and an intramolecular 
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Coulomb repulsion. Since xo can be calculated within 
DMFT, X follows from Eqn. @). 

Putting Ti — T2, T3 ^ Ti and taking the Fourier trans- 
form with respect to ri — T3 in the T — > limit, we obtain 

/•OO 

X'°^(zt^„)= / p'°%e)/{e-iiOn), (4) 
Jo 

where 

p'-(e) = ^ |(n,7V- 2| ^c„,Tc„a|0,iV)P 

x6{e-Eo{N)+E,,{N-2)) + ... (5) 

Here \n, N) is the nth excited state of the system with 
N electrons and the energy £"„ (N) . The term shown de- 
scribes the removal of an electron pair; " ..." indicates the 
addition of an electron pair. 
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FIG. 1. The ratios S/So and x/Xo drop as a function 
of U/W for A = 0.7. For the non-Jahn Teller model A x a, 
these ratios drop rapidly as U increases. In contrast, the pair- 
ing susceptibility for i5 x e is very resistant to increasing U. 
The results were obtained from exact diagonalization for an 
impurity model with five host sites. 

The new physics can be best understood by examin- 
ing a sum rule for the spectral function p'°'^(e). In the 
simplest Jahn- Teller case, E x e, 

/"OO 

/ p^°'{e)de = S^A, (6) 
Jo 

for a half-filled band and in the limit of very large U 
and very small A. In this limit the Jahn- Teller phonons 
produce local singlets on the molecules, 

-ilZ^^T^mil^ac), (7) 
^ ni 

with pairing via the molecular quantum number to. Al- 
though the Jahn- Teller effect competes with intermolec- 
ular hopping, a large U reduces hopping, so that singlets 
form even for small A. In contrast, for pg°'^(e) the corre- 
sponding sum rule gives only S'o = 1. Since x'°'^ tends 



to be larger than Xg""^, the effective interaction r'°'^ (Eq. 
(^)) tends to be attractive. The existence of local sin- 
glets (Eq. (^) means that the probability for removing 
or adding two electrons with the same m quantum num- 
ber is very high. In contrast, for xo the electrons' m- 
quantum numbers are independent and Xo°'^ tends to be 
smaller. 

As U increases, spectral weight is shifted upwards in 
energy, which tends to decrease x'°'^- However, this is 
partly compensated by an increase of the integrated spec- 
tral weight S, since the Jahn- Teller effect wins when hop- 
ping is reduced (see in Fig. |l] ) . Increasing U therefore 
does not rapidly eliminate a negative r'°'^, as one might 
have expected. In contrast, for the A x a model the 
Coulomb repulsion U and the electron-phonon coupling 
directly compete, so S (and therefore r'°'^) drops quickly 
as U is increased. These results illustrate one impor- 
tant aspect of molecular solids with Jahn- Teller phonons: 
counter-intuitively, Coulomb interactions can in certain 
respects actually help electron-phonon coupling. Capone 
et al. Ipzl have reached similar conclusions for A4C60 us- 
ing a different approach. 

Another important aspect is screening. The Coulomb 
interaction is well-screened by the transfer of electrons 
between the molecules Although this helps su- 

perconductivity, it also normally implies an equally ef- 
fective screening of the electron-phonon interaction it- 
self. However, since the Hg phonons do not shift the 
center of gravity of the electronic levels, they cannot 
be efficiently screened by charge transfer on and off the 
molecule @|9|. 

Both these effects are missing for A^ phonons. Ag 
phonons furthermore tend to cause instabilities when 
coupled to a degenerate level. Within a semiclassical 
approximation, a molecular solid with U = Q becomes 
unstable when A > 1/(2A^), where N is the orbital de- 
generacy. A QMC calculation for a T x a model supports 
this result, whereas in the T x h case the system stays 
metallic for A < 1 (and U = Q). To be able to use a 
reasonably large A, we therefore study the A x a system 
below. 

Although it is now clear that Jahn- Teller phonons can 
cause local pairing, as described by x'°'^/Xo°'^ s-'^d r'°'^, 
superconductivity requires the formation of a coherent 
state through the solid. With a finite coherent metallic 
weight at the chemical potential and an attractive inter- 
action F, the divergent unperturbed uniform pair prop- 
agator Xo mediates a pairing instability towards forming 
the coherent superconducting state. 
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FIG. 2. Tc as a function of A according 
to Migdal-Eliashberg (dashed line) and DMFT theories for 
the T X h (o) and Axa (a) couplings at half-filling. The 
parameters are ujph/W = 0.25 and U = 0. The Txh results 
for U = W {•) are also shown. 



To obtain explicit results for the superconducting tran- 
sition, we use DMFT calculations solving the effective 
impurity problem using QMC pO| , ^ . We first discuss 
the case U ~ 0. Fig. || shows Tc as a function of A ac- 
cording to the DMFT and EUashberg H theories. The 
Eliashberg theory is expected to overestimate Tc of doped 
Ceo both because of the violation of Migdal's theorem 
and because the mean-field Eliashberg equations are in- 
sufficient for a bad metal |10|. Surprisingly, in the Txh 
model for U = 0, the Eliashberg Tc remains accurate 
even up to relatively large values of A. In contrast, for 
the Axa problem the Eliashberg theory fails at smaller 
A ||2^. The DMFT calculation shows a maximum Tc at 
A ^ 1 due to a rapid drop of xo with increasing A beyond 
this range (Eq. (|l|)). For small A, xo goes as 1/(1 + A), 
which renormalizes A to A/ (1 -I- A) in the McMillan equa- 
tion p^ ]. For a larger A, however, xo drops much faster in 
the DMFT than in the Eliashberg theory; the formation 
of a coherent state is less efficient, since spectral weight 
rapidly transfers away from the chemical potential as the 
system approaches a metal-insulator transition. 



We next discuss finite U, connecting to Fig. |l| The 
solid points of Fig. || show the overall reduction in Tc for 
finite U. Fig. |^ shows Tc as a function of U for the Txh 
and Axa models. For A x a, Tc drops quickly when U 
increases, as expected. However for T x h, Tc is more 
resistant to increasing U |2^] . This is consistent with the 
local pairing of Fig. ^ and illustrates the importance of 
treating explicitly the dynamic interplay between Jahn- 
Teller phonons and electrons in molecular systems with 



0.02 



T/W 



0.01 




.^Axa (A phonons) 



0.2 0.4 0.6 0.8 1 

u/w 

FIG. 3. Tc as a function of U for the Txh and Axa 
models for half-filling. The parameters are A — 0.6 and 
i^ph/W = 0.25. The figure illustrates the important difi'er- 
ence between Hg and Ag phonons. 
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FIG. 4. Tc as a function of doping n for difi^erent values of 
U for Txh coupling. The parameters are ujph/W = 0.25 and 
A = 0.6. The figure illustrates the strong doping dependence 
for U/W > 0.4. 

Experimentally, Tc drops quickly in fuUerides when 
the doping n is reduced below three electrons per Cgo 
molecule |24|. This cannot be explained within Eliash- 
berg theory: reducing n from 3 induces a slight increase 
of the density of states at the Fermi energy , which 
should increase A and Tc- This has been taken as ev- 
idence for an electron-electron mechanism of supercon- 
ductivity Fig. ^ shows the doping dependence of 
Tc in DMFT. For smaU U, Tc drops slowly until n ^ 2 
or 4 and then starts dropping much faster: r'°'^ drops 
rapidly here, probably because local pairing is inefficient 
once the average number of electrons per molecule drops 
below two. For U/W > 0.4, Tc drops more quickly as 
n = 2 is approached. The system can gain a particular 
large Jahn- Teller energy aX n — 2 |^,^; this moves the 
system towards a metal-insulator transition and shifts 
spectral weight in the one-electron Green's function away 
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from the chemical potential ||T|]. The shift in spectral 
weight rapidly reduces xo and therefore Tc. Thus the 
strong doping dependence can be explained within an 
electron-phonon mechanism, and there is no need to as- 
sume an electronic mechanism. 

For conventional superconductors, retardation is im- 
portant in reducing the effects of the strong Coulomb 
repulsion ||^. For the fuUerides, the retardation effects 
are small, since W ~ u!ph [||J|]- Local pairing and 
screening are therefore crucial in reducing the effects of 
the Coulomb repulsion for the fuUerides. An increasing 
Coulomb interaction does not much damage supercon- 
ductivity, since the concomitant reduction in hopping fa- 
vors the Jahn- Teller pairing. This leads to new physics in 
these strongly correlated low-bandwidth molecular solids. 
The importance of local pairing is consistent with the 
short coherence length, which is only about three times 
the Cgo-Cgo separation p^ , p^ . 
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